Draft version February 5, 2008 

Preprint typeset using I4TgX style emulateapj v. 6/22/04 



in 
o 
o 

(N 

< 

0^ 



> 
00 

o 
in 
o 

:^ 

Oh 

6 



X 



DIFFUSE LIGHT IN THE VIRGO CLUSTER 

J. Christopher MiHos, Paul Harding, John Feldmeier', and Heather Morrison 

Department of Astronomy, Case Western Reserve University, 10900 Euclid Ave, Cleveland, OH 44106 

Draft version February 5, 2008 

ABSTRACT 

We present deep optical imaging of the inner ~ 1 .5°xl .5° of the Virgo cluster to search for diffuse intracluster 
light (ICL). Our image reaches a Icr depth of /iv=28.5 mag/arcsec^ — L5 mag/arcsec^ deeper than previous 
surveys — and reveals an intricate web of diffuse intracluster light. We see several long (>100 kpc) tidal 
streamers, as well as a myriad of smaller-scale tidal tails and bridges between galaxies. The diffuse halo of 
M87 is traced out to nearly 200 kpc, appearing very irregular on these scales, while significant diffuse light 
is also detected around the M84/M86 pair Several galaxies in the core are embedded in common envelopes, 
suggesting they are true physical subgroups. The complex substructure of Virgo's diffuse ICL reflects the 
hierarchical nature of cluster assembly, rather than being the product of smooth accretion around a central 
galaxy. 

Subject headings: galaxies: clusters: individual (Virgo) — galaxies: interactions 



1. INTRODUCTION 

The diffuse intracluster light (ICL) in galaxy clusters is a 
valuable tool for studying their dynamical history. First iden- 
tified (see Vflchez-Gomez 1999 for review) and quantified 
(e.g., Uson, Boughn, & Kuhn 1991; Bernstein et al. 1995; 
Feldmeier et al. 2002, 2004) using deep broadband imaging, 
the ICL has also been detected both in individual stars (Fergu- 
son, Tanvir, & von Hippel 1998; Durrell et al. 2002) and intra- 
cluster planetary nebulae (IPNe; see Feldmeier et al. 2004 and 
Aguerri et al. 2005 and references therein), the latter open- 
ing up a new field studying ICL kinematics (Arnaboldi et al. 
2004; Gerhard et al. 2005). On the theoretical front, N-body 
simulations of galaxy clusters are now starting to yield high- 
resolution predictions of the spatial and kinematic distribu- 
tion of ICL (e.g., Willman et al. 2004; Murante et al. 2004; 
Sommer-Larsen, Romeo, & Portinari 2005), and showing that 
the formation of the ICL is intimately linked to tidal stripping 
during the hierarchical assembly of clusters. The phase-space 
distribution of the ICL thus holds important information about 
the detailed dynamical state and assembly history of a galaxy 
cluster. 

In this context, the Virgo cluster (at an adopted distance 
of 16 Mpc) is of particular interest as it shows many signs 
of being a dynamically complex environment, ideal for the 
production of ICL. The cluster possesses both spatial (Bin- 
gelli, Sandage, & Tammann 1987) and kinematic (Bingelli, 
Popescu, & Tammann 1993) substructure, with different sub- 
groups possessing different morphological mixes of galaxies 
(Bingelli et al. 1987). Many of the cluster galaxies are kine- 
matically disturbed (Rubin, Waterman, & Kenney 1999), in- 
dicative of strong, ongoing tidal interactions. Deep photo- 
graphic imaging (e.g., Malin 1994; Katsiyannis et al. 1998) 
also show tidal features around many Virgo galaxies, as well 
as an extended stellar envelope around M87 (Weil, Bland- 
Hawthorn, & Malin 1997). Studies of the ICL in Virgo also 
hint at an irregular structure: the distribution of intracluster 
PNe show marked field-to-field variations (Durrell et al. 2002; 
Feldmeier et al. 2004; Aguerri et al. 2005), suggestive of a 
poorly mixed ICL. 
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However, quantifying the large-scale structure of the ICL in 
nearby clusters like Virgo has actually proved quite difficult. 
The relatively small field of view of most telescope/CCD sys- 
tems makes it hard to survey nearby clusters of large angular 
size. Deep surface photometry also is hindered by the need 
to accurately flat field over degree scales, while IPNe surveys 
to date have not gone deep enough to uniformly sample the 
full line-of-sight depth of the Virgo cluster Accordingly, the 
distribution of ICL in Virgo is still poorly determined, yet it 
is here that there is a wealth of high spatial resolution ancil- 
lary data on the properties of the galaxies and hot intracluster 
medium. 

To address this disparity and make quantitative, large-scale 
measurements of the diffuse light in Virgo, we have begun 
a deep imaging study of Virgo's core using Case's Burrell 
Schmidt telescope. Because the Schmidt's wide field of view 
is imaged onto a single CCD chip, highly accurate flat fielding 
is possible. Its closed-tube design and extensive optical baf- 
fling significantly reduce scattered light problems compared 
to open-tube telescopes. As such, the Burrell Schmidt is opti- 
mally suited for deep, wide-field surface photometry. Here we 
report on the first season of observations which have revealed 
Virgo's intricate web of diffuse ICL. 

2. IMAGING TECHNIQUE 

Deep images of the Virgo cluster were taken during pho- 
tometric nights in March and April 2004 using the recently 
refurbished 0.6m Burrell Schmidt telescope (for technical de- 
tails of the refurbishment see Harding et al. 2005 in prepara- 
tion, and Feldmeier et al. 2002 for details of the observing and 
data reduction strategies). The field of view of the Schmidt 
and its SITe 2048x4096 CCD is approximately 1.5°x0.75°, 
with 1.45" pixels and the long axis of the chip oriented east- 
west. A total of 72 object images were taken, covering the 
inner 2.25 square degrees of Virgo. Individual images were 
flattened using a night sky flat built from 127 offset sky point- 
ings bracketing the object exposures in time and position. The 
object images and blank skies are all 900s exposures taken in 
Washington M and transformed into Johnson V magnitudes. 
Skv subtraction is a challenge when mosaicking together 
f '^"rJ?8«f'=L^l'pP^CTa?!^^. jtend only a portion of the cluster it- 
self. To remove stars and galaxies from the individual images. 



Fig. 1 . — Diffuse light in the Virgo cluster core. North is up; east is to the left. The white levels saturate at /Uv~26.5, while the faintest features visible have a 
surface brightness of /iv~28.5. 



we first mask stars out to where the PSF falls below 3 ADU 
(Hy=2S), then subtract their extended wings out to 0.5 ADU 
(/xv=30). IRAF's OBJMASK task is then used to mask out all 
extended sources 2.5a above sky, which essentially removes 
all bright galaxies. We then remove sky in a multi-step pro- 
cess. We first rebin the masked images into smaller images 
by calculating the mode in 32x32 pixel bins, and then identify 
regions on the images which are located far away from any 
bright galaxy and considered "pure sky"; on average, about 
5% of each image is covered by pure sky. Then for each im- 
age, we calculate a DC sky level from the mode of the counts 
in these regions and subtract it from the image. This sky 
level ranges from 1 100-1500 ADU (/iv=2 1.6-2 1.3), depend- 
ing on a variety of factors such as hour angle, airmass, time of 
night, and air glow variability. Next, we perform an iterative 
plane-fitting process where the pure sky regions are used to si- 
multaneously fit individual residual sky planes to each binned 
image, constrained to minimize frame-to-frame deviations in 
the pure- sky regions of the final sky- subtracted images. The 
edge-to-edge gradient of the final fitted sky planes is typically 



5 ADU, or ~ 0.35% of the DC sky level.^ 

After sky subtraction of the individual images, the images 
are registered and medianed together to form the final mo- 
saic. At this point one last plane is also fit and subtracted 
from the mosaic using the pure sky regions; the edge-to- 
edge gradient of this final plane is 2 ADU 0.1% of sky). 
We then improve the signal-to-noise by medianing the mo- 
saic in a 9x9 pixel boxes (~ 1 kpc^), which brings out the 
faintest features. In the raw image, the Icr level in the sky 
is 4.5 ADU (/iv=27.5), while in the binned, medianed image 
1(7=1.7 ADU (/iv=28.6). This Hmiting surface brightness is 
approximately 1.5 mag/arcsec^ deeper than the deep photo- 
graphic imaging of MaUn (1994). 

As a final note, one possible source of confusion in our 
study is contamination due to the presence of any galactic cir- 
rus in the field, which will be visible as large-scale diffuse 

^ It is important to recognize that fitting and subtracting sky planes can 
impose a systematic underestimate of the Virgo ICL, since any degree-scale 
diffuse ICI^ component will also be subtracted at this step. 



3 




Fig. 2. — Diffuse light in the Virgo cluster, color coded by surface bright- 
ness: Blue: fiv=25-26; Yellow: /xv=26-27; Green: /Jv=27-28; Red: fiv=28- 
29. North is up; east is to the left. 
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Fig. 3. — Schematic diagram showing location of diffuse features discussed 
in the text. 



reflected light. Examining the IRAS 100//m maps of Virgo, it 
is clear that a ring of cirrus ('^ 1 .5° in diameter) surrounds the 
Virgo core and may contaminate fields to the southwest and 
north of our field. Fortuitously, however, the Virgo core itself 
is quite clear of galactic cirrus, arguing that features we detect 
are bona fide structures in Virgo's diffuse light. 

3. DIFFUSE LIGHT IN VIRGO'S CORE 

A greyscale representation of the final binned mosaic is 
shown in Figure 1, while Figure 2 shows the image colored 
according to surface brightness. The faintest features visible 
are at /iv~ 28.5 or about 0.1% of sky. A wealth of diffuse 
features can be seen in Virgo's ICL, ranging from extended 
low surface brightness envelopes to long, thin streamers, as 
well as smaller-scale tidal features associated with many of 
the Virgo galaxies. Our goal here is to focus on the qualitative 
morphology of the ICL; a more detailed quantitative analysis 
will be addressed in later papers. Figure 3 shows a finding 
chart for the features described below. 

In the image, two long streamers can be seen extending to 
the northwest from M87 at the lower left. One streamer (A) 
projects through the pair of galaxies NGC 4458/61 and ex- 
tends beyond, toward a group of galaxies to the north. There 
is a surface brightness gradient along its length; near M87 it 
has a surface brightness of /iv=27, dropping to /iv=27.8 near 
the pair From the point at which it emerges from M87's stel- 
lar halo to where it fades into the background, the streamer 



has total length of 38' (178 kpc) a characteristic width of 3f5 
(16 kpc), and a total luminosity of ~ lO' L©. Because NGC 
4458 and NGC 4461 have a very high velocity relative to one 
another (Av = 1296 km/s) it is unlikely the streamer comes 
from any strong interaction between the pair; instead, strip- 
ping from one of the galaxies individually is more likely. 

The second streamer (B) is perhaps even more remarkable. 
Also appearing to emanate from M87 to the northwest, it mea- 
sures approximately 28' (130 kpc) long and only 1.2' (6 kpc) 
wide, with a peak surface brightness of /iv=27.0. The low 
surface brightness dwarf VCC 1 149 is projected along the 
stream, just before the position at which the stream blends 
into M87's halo; we may well be seeing the destruction of 
this dwarf in the tidal field of the cluster — a process which 
may contribute to the variation in the faint end slope of the lu- 
minosity function in clusters. The thinness of the stream sup- 
ports this conjecture, as the low velocity dispersion of dwarf 
galaxies ensures that stripped tidal streams will be dynami- 
cally cold, and not diffuse significantly in the tangential di- 
rection. The thinness and linearity of the streamer also argue 
that the dwarf must be on a highly radial orbit, or one that is 
seen projected along its orbital plane. 

Not all detected streams are linear, however A broad, 
curved stream (C) can be seen surrounding a low surface 
brightness dwarf; the curvature of this stream, arcing away 
from M87, suggests that the galaxy is orbiting within or 
around the M86 group (or perhaps the NGC4435/8 pair) 
rather than around M87 itself. 

Aside from long ICL streams, many galaxies show tidal 
features on smaller scales. Perhaps the most obvious is the 
tidal tail to the north of the interacting pair NGC 4435/8 (D). 
Originally identified by Malin (1994) in photographic imag- 
ing, our imaging traces it to fainter surface brightnesses (/iy > 
28), where it takes a sudden 90° bend to the west. This kind of 
tidal "dogleg" is expected during close and slow encounters in 
a galaxy cluster, where the tidal forces of both the individual 
galaxies and the cluster itself conspire to create rather com- 
plex tidal morphologies (Mihos 2005, in preparation). Other 
notable features include NGC 4388's tidal plumes to the west 
and northeast (E), a southern tail in NGC 4425 (F; both noted 
by Malin 1994) which connects to VCC 987, another tail ex- 
tending from IC 3349 (G) and passing through VCC 942, and 
a loop to the east of IC 3475 (H). 

At low surface brightness, the giant elliptical M87 shows 
a very extended and irregular stellar envelope. From pho- 
tographic images, Weil et al. (1997) and Katsiyannis et al. 
(1998) show M87's luminous halo extending to ^ 100 kpc 
scales; on our image it can be traced further out to 175 kpc 
(I). Here the outer isophotes are not simply the radial exten- 
sion of M87's inner regular luminosity profile; traced to the 
northwest, the isophotes become boxier and show irregular 
loops and shells which may be remnants from earlier strip- 
ping from, or cannibalization of, smaller galaxies. The lop- 
sidedness and irregularity of M87's halo is again indicative of 
a diffuse envelope which is not fully relaxed. 

Because of the close proximity of M86, M84, and a number 
of smaller galaxies, uniquely tracing out any individual halo 
in that complex of galaxies proves difficult. The isophotes 
of both M86 and M84 appear fairly regular out to a surface 
brightness of /iy-- 26.5 (J), as noted by Malin (1994). At 
fainter surface brightnesses, the isophotes blend together and 
become quite irregular, and we see a "crown" of diffuse light 
(at [iy= 27-28) to the north of the M84/M86 system (K), and a 
diffuse extension to the south of the pair, blending into the net- 
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work of galaxies to the south. Whether this extended envelope 
is truly physical or simply a projection effect remains unclear. 
Also visible around M86 are several low surface brightness 
filaments (L1-L3), one to the northeast (/iv~ 25.5), one to 
the west (/Uv~ 26.5), and one to the north (/Uv~ 27.5). 

We also see a common envelope of light (M) surround- 
ing the complex of galaxies to the south of M86, including 
NGC 4413, IC 3363, and IC 3349, along with a tidal bridge 
connecting NGC 4413 and IC 3363, and perhaps extending 
to NGC 4388 to the east. This envelope has fairly irregular 
isophotes, arguing it is more than just simply the extended 
profiles of galaxies overlapping in projection. This group of 
galaxies may be a physically coherent subgroup within the 
Virgo cluster. 

Finally, Oosterloo & van Gorkom (2005) have recently 
identified an extended (110x25 kpc) HI cloud streaming from 
NGC 4388 which they attribute to ram-pressure stripping due 
to the lack of any stellar counterpart in previous photographic 
imaging. We also find no extended counterpart on our im- 
age down to a surface brightness limit of /xv=28.5, save for 
faint streamers close to the main body of NGC 4388 coinci- 
dent with the extended emission Une region noted by Yoshida 
et al. (2002). 

4. DISCUSSION 

Our deep imaging of the Virgo cluster has revealed a com- 
plex network of extended tidal features, indicative of the on- 
going stripping and tidal evolution of galaxies in Virgo. The 
ICL is not radially symmetric around M87, the nominal cen- 
tral galaxy of Virgo; much of the diffuse light is centered upon 
the M84/M86 complex. The impact of hierarchical assembly 
of clusters on the ICL is clear here — rather than ICL growing 
simply via smooth accretion around a central galaxy, it distri- 
bution reflects the substructure inherent in the cluster, as ev- 
idenced in cosmologically-motivated N-body simulations of 
ICL (e.g., Dubinski 1998; Willman et al. 2004, Rudick et al. 
in preparation). Characterizing the ICL in terms of a simple 
radial profile will prove misleading in all but the most regular, 
relaxed clusters. 

It is also instructive to compare our map of the ICL in 
Virgo's core to estimates of the luminosity density derived 
from studies of planetary nebulae by Feldmeier et al. (2004) 



and Aguerri et al. (2005). Because IPNe detection depends 
strongly on the underlying surface brightness, IPNe studies 
have complicated selection and completeness functions, par- 
ticularly in fields which cover both bright galaxies and diffuse 
intracluster space. With this in mind, we will restrict the dis- 
cussion here to a quaUtative comparison, and leave a more 
quantitative analysis to a future paper 

Both IPNe studies show that the IPNe density is highly vari- 
able in the Virgo core; our image confirms this complexity 
and variation in the underlying ICL. The highest IPNe densi- 
ties are found in Feldmeier et al.'s "Field 3," which our image 
shows is completely embedded in M87's diffuse halo. Much 
of the IPNe detected in this field are likely to be associated 
with M87 itself; radial velocity measurements of the IPNe in 
this field agree with this assessment (Arnaboldi et al. 2005). 
Aguerri et al.'s "SUB" field also shows a relatively high PNe 
density; aside from containing the luminous ellipticals M84 
and M86, this field also partly covers the diffuse envelope 
surrounding the group of galaxies to the south of M86. The 
lowest IPNe densities are reported in Aguerri et al.'s "LPC" 
field, which also corresponds to the lowest luminosity density 
in our image; indeed, over most of this field, we detect no 
light above the background. 

The qualitative match between the luminosity density 
traced by IPNe and broadband hght argues that IPne are effec- 
tive tracers of the ICL in galaxy clusters, and vice- versa. Our 
deep imaging provides a "finding chart" for IPNe searches; 
followup spectroscopy can then be used to study the kine- 
matics of discrete structures such as tidal tails and diffuse 
envelopes around galaxies. Our survey of diffuse light in 
Virgo continues: we have follow-up imaging of additional 
portions of the Virgo cluster, and subsequent papers will focus 
on quantitative descriptions of Virgo's ICL, as well as a sur- 
vey of extremely low surface brightness galaxies in the cluster 
core. 
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